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Abstract 

The remnant G347.3-0.5 exhibits strong shell emission in the radio and X-ray bands, 
and has a purported detection in the TeV gamma-ray band by the CANGAROO-II 
telescope. The CANGAROO results were touted as evidence for the production of 
cosmic ray ions, a claim that has proven controversial due to constraining fluxes as¬ 
sociated with a proximate unidentified EGRET source 3EG J1714-3857. HESS has 
now seen this source in the TeV band. The complex environment of the remnant 
renders modeling of its broadband spectrum sensitive to assumptions concerning 
the nature and parameters of the circumremnant medium. This paper explores a 
sampling of reasonable possibilities for multiwavelength spectral predictions from 
this source, using a non-linear model of diffusive particle acceleration at the shocked 
shell. The magnetic field strength, shell size and degree of particle cross-field diffu¬ 
sion act as variables to which the radio to X-ray to gamma-ray signal is sensitive. 
The modeling of the extant data constrains these variables, and the potential impact 
of the recent HESS detection on such parameters is addressed. Putative pion decay 
signals in hard gamma-rays resulting from hadronic interactions in dense molecular 
clouds are briefly discussed; the requisite suppression of the GeV component needed 
to accommodate the 3EG J1714-3857 EGRET data provides potential bounds on 
the diffusive distance from the shell to the proximate clouds. 
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1 Introduction 


Supernova remnants (SNRs) are believed to be the principal candidates for 
generating galactic cosmic rays, relativistic ions and electrons with energies 
from the GeV range up to 10^^ eV and beyond. A Holy Grail of cosmic ray 
physics has been to secure a radiative detection that would unambiguously 
indicate the presence of cosmic ray ions in proximity to one or more supernova 
remnants, thereby establishing a direct connection between source and diffuse 
population. The search has therefore focused on gamma-ray emission signaling 
the decay of neutral pions, 7r° ^ 77; these are presumably generated in 
collisions between cosmic ray protons accelerated at the shocked outer shells 
of remnants and ambient hydrogen. This is a difficult task, due largely to the 
relatively low fluxes expected, which approach the limits of sensitivity in the 
100 MeV band for past hard gamma-ray experiments such as EGRET. 

While this will be addressed by the upcoming GLAST mission, scheduled 
for launch in 2007, the Imaging Atmospheric Cerenkov Telescope (lAGT) tech¬ 
nique has focused a large amount of time and energy over the last decade in 
searching for pion decay signatures in the TeV band. Seminal steps along this 
path were the announcement of detections of the young remnant SN1006 by 
GANGAROO (Tanimori et ah 1998), and even younger Gas A by HEGRA 
(Aharonian et al. 2001) in the northern hemisphere. A second southern rem¬ 
nant, G347.3-0.5 (radio designation, or RX J1713.7-3946 in X-rays), has a 
claimed detection by the GANGAROO team (Enomoto et al. 2002) in its 
northwest rim, an asymmetry mirroring the GANGAROO TeV observation 
of SN1006 (which contrasts its beautifully symmetric radio and X-ray maps: 
e.g., see Koyama et al. 1995). The RX J1713.7-3946 discovery paper spawned 
a mini-controversy over whether or not the detection could be interpreted as 
evidence for hadrons in the remnant. The model interpretation of this remnant 
forms the focus of this work. 

The GANGAROO team offered a hadronic interpretation for the signal 
they claimed in Enomoto et al. (2002), using particle distributions that would 
result from linear shock acceleration models. This assessment was disputed 
in subsequent papers in the literature, in particular Reimer & Pohl (2002) 
and Butt et al. (2002). These dissenting opinions argued that in order to fit 
the spectrum published in Enomoto et al. (2002) in the TeV range, the model 
pion decay flux at 100 MeV energies would be so large that it would far exceed 
that of the EGRET unidentified source 3EG J1714-3857 lying in the field of 
view of G347.3-0.5. Whether or not the EGRET source is associated with the 
remnant is immaterial: it can serve as a putative upper bound, suggesting that 
an inverse Gompton model using cosmic microwave background (GMB) seed 
photons would be a more viable explanation. 
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This paper offers additional perspectives on this source, specifically explor¬ 
ing predictions of radiation models that incorporate the physics of non-linear 
acceleration at remnant shocks. Essentially it updates the careful exploration 
of Ellison, Slane and Gaensler (2001) that predated the full publication of the 
CANGAROO data on this source, providing elements of a parameter survey. 
An update on radio and X-ray observations for this remnant has recently been 
provided by Lazendic et al. (2004). Yet crucial to multi-wavelength modeling 
is a recent observational development: the new stereoscopic HESS telescope 
in Namibia has now detected G347.3-0.5 with an image that nicely matches 
the X-ray map, albeit with a slightly lower flux and a flatter spectrum from 
the GANGAROO claim. The HESS results are presented in Aharonian et al. 
(2004), and complicate the modeling picture considerably, particularly since 
HESS has not detected SN 1006 (Aharonian et al. 2005), providing upper lim¬ 
its an order of magnitude or more below the GANGAROO data. The model 
results presented here for G347.3-0.5 will leave the validity of the different 
TeV data sets to further appraisal by the observational community. 


2 Multi-wavelength Models from Non-linear Shocks 


Model predictions for supernova remnant emission often invoke linear mod¬ 
els of acceleration, where the accelerated particles have no influence on the 
hydrodynamics of the expanding shock (e.g. Drury, Aharonian & Volk 1994; 
Gaisser, Protheroe & Stanev 1998; Stumer et al. 1997; Enomoto, et al. 2002; 
Reimer & Pohl 2002). Such linear invocations generally suffice to predict X-ray 
and TeV fluxes to within factors of 3 — 10 . Usually, acceleration is considered 
to take place at the expanding forward shock, though increasingly the impor¬ 
tance of acceleration at the reverse shock is being realized (e.g. Decourchelle, 
Ellison & Ballet 2000; Ellison, Decourchelle & Ballet 2005). Acceleration in 
strong forward shocks can be very efficient just prior to, or in the early stages 
of, the Sedov phase, leading to a sizeable fraction of the total energy bud¬ 
get being placed into cosmic rays. Accordingly, energy flux conservation in 
the shock hydrodynamics is profoundly influenced by the contribution from 
the highest energy cosmic rays, modifying the shock compression, and hence 
inducing a non-linear feedback on the acceleration process. 

Such non-linear acceleration effects are well-documented in the cosmic ray 
literature (e.g. see Drury 1983; Jones & Ellison 1991 for reviews). Their key 
characteristic, a direct consequence of diffusive spatial scales generally be¬ 
ing an increasing function of particle momenta, is a deviation from pure 
power-law behavior in the ion and electron distributions, with an upward 
curvature indicating more efficient acceleration of the most energetic parti¬ 
cles. Such modifications to simpler, linear (test particle) models can influ¬ 
ence radiative predictions for comparison with data by factors of 3 — 10 or 
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more, and therefore become imperative for the theory/observation interface. 
Detailed models of non-linear acceleration come in different varieties: Monte 
Carlo simulations have been documented in Ellison & Eichler (1984), Jones 
& Ellison (1991) and Ellison, Baring & Jones (1996), and analytic/numerical 
techniques involving diffusion-convection equations have been developed by 
Kang & Jones (1995), Berezhko, Yelshin & Ksenofontov (1996), Maikov (1997) 
and Blasi (2002). Extensively-developed applications of these complementary 
approaches to SNR emission yield similar radiative predictions for similar pa¬ 
rameters, so that there is a robustness in the theoretical modeling (e.g. see 
Baring 2000, for a discussion). 

In this paper, a rather compact analytic representation of the non-linear 
problem developed for ions in Berezhko & Ellison (1999), and extended to 
treat electrons in Ellison, Berezhko & Baring (2000), is used. This formulation 
models the non-linear feedback and uses piecewise continuous broken power- 
laws to approximate the non-thermal distributions. The indices of the power- 
law segments are determined self-consistently, and the resulting distributions 
are extremely facile for folding into radiation models for remnants. The reader 
is referred to Ellison, Berezhko & Baring (2000) for a comprehensive discussion 
of the acceleration theory, and Ellison, Slane & Gaensler (2001) for a compact 
presentation on its application to G347.3-0.5. 


2.1 Leptonic Spectral Models 


In the light of the gamma-ray constraints imposed by the EGRET source in 
the locale of G347.3-0.5 on the sky, it is prudent to focus on inverse Gompton 
models for TeV photon production, following Reimer & Pohl (2002), and many 
authors in applications to other remnants. The radio to X-ray continuum is 
modeled with synchrotron emission from a shock-accelerated distribution of 
electrons, and the same distribution upscatters GMB photons to gamma-ray 
energies. In such models, spectroscopically, ions go along for the ride, their 
major role being a power drain for radiative efficiency. Photon production 
rates are standard calculations, assuming isotropy of electrons, and the reader 
is referred to Baring et ah (1999), for example, for details. 

Ellison, Slane & Gaensler (2001) used the hydrodynamics of a wind-shell 
interaction to dehne shock parameters, using these as input into the non-linear 
acceleration formalism of Ellison, Berezhko & Baring (2000) that was then 
used to generate photon emission models. Here, the hydrodynamic consider¬ 
ations were not explicitly made, but parameters similar to those of Ellison, 
Slane & Gaensler (2001; specihcally a hybrid of their models A, B and G) 
were adopted to illustrate the main features. Spectra from the hrst sample 
model are illustrated in Fig. 1, together with relevant, published data. This 
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is a model with a modest upstream field of = 2/i G and compression ratio 
Got = 5. The field compression of course depends on the shock obliquity, so 
here we essentially model an almost perpendicular shock, so that = 10/i G. 



Logio[Energy, (MeV)] 

Fig. 1. The multi-wavelength spectra for synchrotron (solid curve) and inverse 
Compton (dashed) components resulting from a non-linear shock acceleration model 
of SNR G347.3-0.5 (RX J1713.7-3946). Model/source parameters are as labeled, 
and are discussed in the text. The observational data consists of 1.36 GHz radio 
observations by ATCA, a 2-10 keV ASCA X-ray spectrum of the northwest rim 
(both detailed in Ellison, Slane & Gaensler 2001), a band at > 300 GeV depicting 
the steep CANGAROO-II spectrum announced in Enomoto et al. (2002), and a band 
representing the EGRET unidentified source 3EG J1714-3857, as taken from Reimer 
& Pohl (2002). Low model densities, Up = 0.02cm“^, render the bremsstrahlung 
and pion decay contributions insignificant at X-ray and TeV energies, respectively. 

The upstream particle density is Up = 0.02 cm“^ for this example, low 
enough that a pion decay signal would be dominated by the inverse Gomp- 
ton one in the TeV band (see Fig. 3 of Ellison, Slane & Gaensler 2001). For 
the same reason, the bremsstrahlung contribution would he far below that 
of synchrotron emission in the X-ray band, and so is also not exhibited. A 
source distance of 6 kpc and a remnant age of fsNR = lO^yrs were adopted. 
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following the choice of Ellison, Slane & Gaensler (2001) that were based on 
the contentions of Slane et ah (1999). Debate over the actual distance to SNR 
G347.3-0.5 will be addressed below. 


The spectra in Fig. 1 differ slightly from those exhibited in Fig. 3 of Ellison, 
Slane & Gaensler (2001), most notably in being a little sharper in the severity 
of the X-ray and TeV turnovers, and in the prominence of the spectral breaks 
in the infra-red and the 100 MeV band. Here, the turnovers are structured, 
broken power-laws, as opposed to the exponential ones opted for in Ellison, 
Slane & Gaensler (2001). Also, the spectral breaks image the structure of 
the underlying electron distributions from the simplihed, analytic, non-linear 
acceleration model. Apart from these subtleties, the results are similar, and it 
is clear that the model predicts a TeV turnover too high to accommodate the 
steep GANGAROO spectrum. 


The global structure of the X-ray and TeV gamma-ray spectral peaks can 
be quickly understood to be dependent simply on the shock speed, and the 
strength of the magnetic field. For the parameter regimes chosen to model SNR 
G347.3-0.5, sufficiently strong B is invoked for the acceleration of electrons 
(not ions) to be limited by synchrotron cooling. In such cases, the maximum 
energy of electrons for remnants in the Sedov phase is (for R > 3/i G, e.g., see 
Baring et ah 1999): 
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Here rj = X/vg is the ratio of the diffusive mean free path to a particle’s gyro- 
radius, a parameter that couples to diffusive transport across helds and that 
determines the acceleration rate of electrons and ions. This form guarantees 
that the X-ray synchrotron turnover energy (i.e., just above the X-ray spectral 
peak in the uFi, representation), which is proportional to , is inde¬ 

pendent of the field strength, and depends only on the shock speed (e.g. see 
de Jager et al. 1996, and the discussion in de Jager & Baring 1997) through 
its dependence on Up and tsNR • 

In contrast, the inverse Gompton (IG) peak occurs at an energy oc R^axe 
that does scale with the held Bd '■ higher helds are needed to move the turnover 
to sub-TeV energies. At the same time, the normalization of the X-ray and 
TeV peaks is controlled purely by the relative energy densities in the soft GMB 
photons and the downstream magnetic held, an elementary characteristic (e.g. 
see Rybicki & Lightman 1979), since the same electrons are emitting in each 
band. Higher helds enhance the synchrotron component relative to the IG one, 
i.e. suppress the gamma-rays if the volume hlling factor / (here it is around 
/ = 0.3) is chosen to hx the radio to X-ray continuum. Hence, it is quickly 
established that, in the representation of Fig. 1, where the y-axis scales 
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with the emitted power, the X-ray and TeV gamma-ray peaks will be of the 
same height when = 3.3/iG. Adjustments in , f and Up would be 
needed to try to match the CANGAROO data in Fig. 1. The very recent 
HESS detection provides additional challenges, as discussed below. 


2.2 A Model for Nearby Distances 


A crucial ingredient in the model presented in Fig. 1 is the assumption that 
the source is 6 kpc distant. Koyama et al. (1997) obtained a distance estimate 
oi D = I kpc based on the column density inferred from soft X-ray absorption 
in ASGA data, when compared to the Galactic Genter direction absorption 
characteristics. This contention has recently been supported by XMM obser¬ 
vations (Gassam-Ghenai et al. 2004) and also evidence for molecular cloud 
material at 1 kpc (Fukui et al. 2003). Glearly, for the observed expansion rate, 
the age of the remnant is therefore uncertain by a factor of 6 (though it would 
still be inferred to be in the Sedov phase), and so also is the inferred mean 
density in the remnant, because of the relationship between density Up and 
distance D for hxed source synchrotron flux: F oc UpD = constant. Slane et 
al. (1999) favored the larger 6 kpc value, principally because they observed 
that GO intensity mapping suggested a projected association of the remnant 
with three molecular clouds that appeared to have GO and HI column densi¬ 
ties matching the X-ray values. The radial velocities of these clouds indicates 
a minimum distance of D = 6 kpc when folded into a galactic rotation model. 
This reasoning obviously depends on the validity of the association. 


The debate in the literature over this distance has not been conclusively 
resolved. Accordingly, presenting a model for a D = 1 kpc case is clearly 
motivated. This translates to assuming a density around Up ~ 0.2 cm“^, for 
which we will not treat bremsstrahlung or line emission here. From Eq. (1), 
this change in density would force i7max,e lower so as to move the synchrotron 
peak below the X-ray band. Adjusting Bd upwards is to no avail, since the syn¬ 
chrotron peak energy is independent of held strength in this cooling-dominated 
domain. Hence to ht the ASGA continuum, lower helds are desirable to emerge 
from cooling-limited acceleration phase space, and thereby render Emax.e sen¬ 
sitive to the choice of Bd ■ For helds R ^ 3/r G, synchrotron cooling is generally 
insignihcant, and the maximum energy of accelerated electrons is controlled 
by the spatial scale of the remnant (e.g. Baring et al. 1999): 
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For the higher density and younger age dictated by choosing D = 1 kpc, a 
held of Bd = 2jj,G yields a synchrotron z/F7 peak in the classical X-ray band. 
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A model spectrum corresponding to this case is depicted in Fig. 2. The 
compression ratio of rtot = 6.4 was slightly higher than in the first exam¬ 
ple (due in part to a higher electron-to-proton abundance ratio and a higher 
value of Fornax,e that spawns increased losses in the non-linear acceleration 
model), and again represented Bd/Bu as in a quasi-perpendicular shock. The 
maximum electron (and proton) energy was 60 TeV in this case, markedly 
higher than in Fig. 1. Consequently, the inverse Compton z/Fh peak moves to 
above 10 TeV, a gross mismatch to the putative CANGAROO points. The flux 
in the IC peak now exceeds that for the X-ray synchrotron peak, since now 
the magnetic held energy density is inferior to that of the cosmic microwave 
background. 



Logio[Energy. (MeV)] 

Fig. 2. The multi-wavelength model spectrum resulting from a non-linear shock 
acceleration model of SNR G347.3-0.5 (RX J1713.7-3946) for an assumed distance 
of Zf = Ikpc, following the proposition of Koyama et al. (1997). The two spectral 
components are synchrotron emission (solid curve) and inverse Compton scattering 
(dashed). Model and source parameters are as labeled. The observational data are as 
in Fig. 1, with the addition of the spectrum for the TeV detection recently 

announced by the HESS telescope team (Aharonian et al. 2004); see text for a 
discussion on normalization. 




Also marked on the Figure is the spectrum announced by the HESS team in 
Aharonian et al. (2004). This spectrum is very different from the CANGAROO 
one. Not only is it flatter ( dn/dE oc as opposed to dn/dE oc E ~‘^-^), it 

possesses slightly higher overall flux. The HESS result is integrated over the 
entire remnant, and therefore differs in normalization from the CANGAROO 
result, which was for the northwest rim portion only. For the purposes of the 
model comparison, the HESS spectrum from just this conhned portion of the 
rim is somewhat lower than the CANGAROO data, reflecting the sensitive 
nature of the HESS telescope array. It is noteworthy that the HESS detection 
is spatially resolved, mirroring the ASCA sky map more closely due to its 
improved angular resolution. Clearly the HESS results are at odds with the 
CANGAROO data, a topical issue that will be the subject of discussion by 
the TeV astronomy community in the coming years. Moreover, the HESS data 
are just as difficult to reconcile with the model spectrum, as depicted. 

This multiwavelength modeling could be viewed as another argument against 
adopting D = Ikpc for the distance to SNR G347.3-0.5. However, it should 
be noted that the models were restricted to the Bohm diffusion limit rj = 
X/vg = 1, where a particle’s diffusive scale is the smallest possible, i.e. its 
gyroradius, corresponding to strong field turbulence. Therefore, less turbulent 
shock environments can increase the diffusive scale, thereby reducing the ac¬ 
celeration rate and lowering Emax,e (see Eqs. (1) and (2) ) without necessarily 
altering other environmental parameters such as Bd and Up. Modihcations 
along these lines can, with a slight adjustment upwards in held strength, ap¬ 
proximately match the CANGAROO data in Fig. 2, but might have greater 
difficulty in accommodating the portion of the HESS signal coming from the 
northwest rim, which would require higher helds in order to reduce the inverse 
Compton to synchrotron ratio at their respective peaks. 


2.3 Hadronic Models? 


Having explored the case for purely leptonic models, it is worth noting that 
hadronic models for the generation of gamma-rays are in principal tenable for 
SNR G347.3-0.5, though not in the simplest form envisaged by Enomoto et 
al. (2002) in their original claim. The objections of Reimer & Pohl (2002) 
and Butt et al. (2002) still stand. However the complex environment of the 
remnant admits the possibility that pions can be created using molecular cloud 
nucleons as targets for the energetic cosmic rays, and these would be somewhat 
remote from the expanding outer shell of SNR G347.3-0.5. Such a departure 
from a one-zone scenario is quite reasonable in these circumstances. 

A spatial separation of accelerator and target alters the expected spectrum 
for pion decay emission. Since the diffusive scale of the higher energy ions is 
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larger than that for the low energy ones that contribute to the decay peak 
at 67 MeV (e.g. see Jones & Ellison 1991; Baring et al. 1999; Berezhko & 
Ellison 1999, for a discussion of such diffusion lengths), the latter ions have 
great difficulty in convecting far upstream of the shocked shell. Consequently, 
the pion signal would be strongly suppressed near 67 MeV in this scenario, as 
would all such emission from particles whose energies correspond to diffusive 
scales less than the distance between the shell and the target cloud. This would 
skew the 7r° decay spectrum to render it peaked near the TeV band in the 
Z/T7 representation, thereby evading the EGRET source bounds. 

Clearly, for this geometrical separation to have a significant impact on the 
pion decay continuum, the distance between remnant shell and target molecu¬ 
lar cloud must be at least a sizable fraction of the radius of the remnant. Such a 
model is not unduly contrived, but does require additional parameters relating 
to spatial scales and turbulence (i.e. rj = X/vg) that are poorly constrained by 
observational data, particularly given the impact of line-of-sight ambiguities to 
the environmental geometry. Furthermore, these parameters might well have 
to vary signihcantly across the remnant in order to reconcile the molecular 
cloud sky geometries with the X-ray and TeV emission morphologies. Finally, 
it might prove difficult to distinguish this pion decay scenario spectroscopi¬ 
cally from inverse Compton emission, given that both would be extremely flat 
and statistics-limited somewhat below the TeV peak. 


3 Conclusions 


In closing, two multi-wavelength models are presented here for broadband 
emission in SNR G347.3-0.5, using fundamental characteristics of non-linear 
shock acceleration as their basis. This offering is not designed to hnd the “per¬ 
fect” model £t, but rather to represent samples from an array of possibilities 
that indicate how such multi-wavelength approaches can constrain model pa¬ 
rameter space. While the examples are focused on an inverse Compton origin 
of TeV emission, which is the preferred interpretation in a low density rem¬ 
nant environment, it is also emphasized that hadronic models with pion decay 
emission in the TeV band are not presently ruled out. Using nucleonic targets 
in proximate molecular clouds that are distinct from the remnant shell can 
potentially generate acceptable model fits to TeV datasets. 

The study of this remnant is presently burdened with the controversy 
concerning the TeV data, an issue also for SN 1006, which has so far not been 
detected by HESS (Aharonian et al. 2005). The contradictory results for SNR 
G347.3-0.5 between CANGAROO and HESS obviously need to be resolved 
in order to propel theoretical understanding at a more rapid pace. This is a 
challenge that the TeV community met in the case of the Crab nebula, and 
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northern hemisphere blazars, that established the branch of TeV astronomy 
on a hrm basis. For supernova remnants, this bridge has yet to be crossed. 
Anticipating watershed results from the increasing number of lACTs to be in 
operation in both hemispheres over the next few years, such a development 
will hopefully be around the corner. 
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